A c=— 2 boundary changing operator for the Abelian sandpile model 
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We consider the unoriented two-dimensional Abelian sandpile model on the half-plane with open 
and closed boundary conditions. We show that the operator effecting the change from closed to 
open, or from open to closed, is a boundary primary field of weight —1/8, belonging to a c = —2 
logarithmic conformal field theory. 
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INTRODUCTION 

The sandpile models are among the simplest and most 
studied non-equilibrium models showing criticality. It 
was introduced by Bak, Tang and Wiesenfeld Q as a 
prototype of a class of dynamical models with generic 
criticality. 

The model we consider here is the standard, unoriented 
two-dimensional Abelian sandpile model (ASM), which 
we first briefly recall. We refer to the recent reviews |g, ||] 
for a more complete account. The model is defined on a 
L x M square lattice. At each site i, there is a random 
variable hi , which counts the number of grains of sand at 
i, A stable configuration is a set of hi € {1,2,3,4}. 

The discrete time evolution is defined as follows. First, 
one grain of sand is dropped on a random site of the 
current configuration Ct, producing a new configuration 
C' t which may not be stable. If C' t is stable, we simply 
set C t +i — C' t . If C[ is not stable (the new hi is equal to 
5), C' t relaxes to a stable configuration Ct+i by letting all 
unstable sites topple: a site with height hi > 5 loses 4 
grains of sand, of which each of its neighbours receives 
1 , something we can write in the form hj — > hj — A^- for 
all sites j. A is called the toppling matrix, equal to the 
discrete Laplacian. Relaxation stops when no unstable 
site remains; the corresponding stable configuration is 
Ct+i- 

Two natural boundary conditions may be imposed 
along the boundaries: open and closed. A boundary site 
is open if it loses 4 grains in a toppling, his three (or 
two) neighbours receiving each 1 grain. It is closed if the 
number of grains it loses is exactly equal to the number 
of neighbours. Thus an open site is dissipative, whereas a 
closed site is conservative (like all bulk sites). In all cases, 
the toppling updating rule hj — > hj — A^- applies, with 
A appropriately defined. The dynamics is well-defined 
provided that not all sites are conservative. 

Dhar gave in Q a first detailed analytical treatment 
of the model. Under very mild assumptions, he showed 
that on a finite lattice, there is a unique probability mea- 
sure P* on the set of stable configurations, that is invari- 
ant under the dynamics. Moreover, P* is uniform on 



its support, formed by the so-called recurrent configura- 
tions. The partition function, defined as the number of 
recurrent configurations, is equal to det A, and goes like 
(3.21) iM . 

It is widely believed that at least some of the critical 
properties can be accounted for by a conformal field the- 
ory (CFT), despite the fact that the ASM show intrinsic 
non local features. The CFT picture has been assumed 
and used by various authors, but very few detailed, ex- 
plicit comparisons have been made. 

Our aim in this Letter is two-fold. First, we add strong 
support to the view that CFT is an appropriate descrip- 
tion by comparing boundary CFT predictions with ASM 
calculations. Second, and because the relevant CFT is 
logarithmic with central charge c = — 2, the ASM pro- 
vides a concrete, lattice realization of such a non-unitary 
theory. In this respect, the cylinder partition functions 
computed below may shed some light on the structure 
of boundary states in logarithmic conformal theories, a 
topic of current interest || . 



THE CONFORMAL FIELD THEORY 

Various results have pointed to a relationship between 
the 2d ASM and the q — > critical Potts model || or the 
c = — 2 explicit Lagrangian realization J2], 
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where 6 and 8 are scalar Grassmanian fields. The rel- 
evance of this field theory was confirmed in Q where 
the unit height variables and other local cluster variables 
were indeed shown to go in the scaling limit to com- 
binations of #'s and derivatives (the off-critical regime 
was also investigated in H], where a relationship with 
the massive extension of (||) was established). 

From the CFT point of view, c = —2 is special as it 
is the simplest example of a logarithmic CFT (and is 
necessarily non unitary) ||. It has been discussed by 
many authors, from different points of view, and gives 
rise to many subtleties. For that matter, we refer the 
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reader to the recent Tehran lecture notes by Flohr jjj 
and Gaberdiel [jllj, and to the references therein. 

The free theory (|l]) has a W-algebra generated by three 
dimension 3 fields and underlies a c = — 2 conformal the- 
ory, which is rational with respect to the extended alge- 
bra Jll], Q. It contains six representations: V_i/g, Vq, 
V3/8 and Vi are irreducible representations constructed 
out from ground states with conformal weight given by 
the subscript, and two reducible but indecomposable rep- 
resentations TZq and TZi, each containing two ground 
states with zero conformal weight. All together they 
form a closed chiral fusion ring, but the characters of 
V_x/8) V3/8, TZq and IZi only are closed under modular 
transformations (and form a fusion subring). The chiral 
characters are computed in 13 TJ] 
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where the 0i are the standard Jacobi theta functions, and 
77 is the Dedekind function fl5[| . 

PARTITION FUNCTIONS 

We want to compute Z2 orbifold partition functions 
for the free theory (|l|) on a cylinder, of perimeter L and 
length M, and modulus t = ijr?. Along the closed 
loop, we allow the two fields 0, to be either periodic or 
antiperiodic, and then sum over the two monodromies. 
On each boundary, we choose either open (Dirichlet) or 
closed (Neumann) boundary condition. 

The calculations are fairly straightforward since the 
functional integral is just the determinant of the Lapla- 
cian, subjected to given monodromy and boundary con- 
ditions. For mixed boundary conditions, a simple calcu- 
lation yields 



det-A 



P/A. mixed 



n 



4M 2 



1,(5) 



with e = 0, i for periodic resp. antiperiodic monodromy. 
These infinite products have been computed in J16[ 



det — Ap, 



ed = 
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mixed 
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The orbifold partition function, defined as half the sum 
of the two determinants, reduces to a single character of 
the triplet algebra at c = — 2 
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When L goes to infinity (q — > 0), the cylinder becomes 
an infinitely long strip with open and closed boundary 



condition on either side. The strip can be conformally 
mapped onto the upper half-plane, with the open and 
closed boundary conditions on either side of the origin, 
on the real axis. The change of boundary condition at 
the origin can then be seen as resulting from the action 
of the ground state of V_i/g |P7[ |. Thus the operator 
0op,ci _ ^cl,op sw it cn ing a boundary condition between 
open and closed is a primary field of weight — |. In the 
realization ([!]), this field is non local in 0, 0. 

Similarly one finds the partition functions for identical 
boundary conditions at both ends of the cylinder 



J open, open 



Xv (?) . 
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Vo has the identity f2 as unique ground state, but 7£o D 
Vo has two, namely Q and its logarithmic partner lu, pro- 
portional to :00:. 

The above partition functions represent the finite size 
corrections to the large volume lattice partition func- 
tions, in the limit L, M — > oo with jj — > — 2ir, i.e. after 
the diverging bulk and boundary free energies have been 
subtracted. For what follows, it will be useful to know 
the leading term of the boundary free energies. For this 
purpose, it is more convenient to consider the partition 
function on a rectangle L x M (i.e. with no periodicity), 
with two edges open and two closed (the sides of length 
M are closed say) . The partition function is the determi- 
nant of the lattice Laplacian with appropriate boundary 
conditions (it can be seen as a lattice version of ([!]), or as 
the actual partition function of the 2d ASM). One finds 



L-l M 



m— n=l 

The Euler-McLaurin formula can be used to evaluate the 
logarithm, and yields (see j^] for similar calculations) 

4G r- 
log Z, 2„ P c„ . = L(M + 1) (L + M + l) log(l + s/2) 

V 2 closed 1 

+ilogM + log23 7 ? (2r) + ... (10) 

where the dots stand for terms that vanish in the limit 
L, M — > oo, and where G = 0.915965... is the Catalan 
constant. 

The coefficient — = log 3.21... of the most diverg- 
ing term is the bulk free energy density, whereas the 
terms linear in L and M give the boundary free energy 
densities (remembering that the number of bulk sites is 
(£-2)(M-2)) 

6G 

/open = — - \ l0g(l + \/2) , 
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Therefore an open site has an excess of free energy equal 
to — with respect to a closed site. The entropy per site 
(= e^) is then equal to 3.21 for a bulk site, 3.70 for an 
open boundary site and 2.07 for a closed boundary site. 
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THE 2D ASM ON THE UPPER HALF-PLANE 

We consider the ASM on the upper half-plane, x S 
Z, y > 1. Depending on the boundary condition, open 
or closed, the ASM is defined in terms of an infinite top- 
pling matrix A op or A c i. Away from the boundary, they 
coincide with the Laplacian on the plane; at a boundary 
site, their action is 

(A op - 4)/(x, 1) = (A cl - 3)/(x, 1) 

= -f(x, 2) - f{x - 1, 1) - f(x + 1, 1). (12) 

We will make an extensive use of their inverse, G op 
and G c i. By the method of images, one easily sees, for 
i = (xi,yi), j = (afejjfe), that (x = x 2 - Xi) 

Go P (i;j) = G(x,y 2 -yi)-G(x,y 2 +y 1 ), (13) 
G c i(i;j) = G(x,y 2 -y 1 ) + G(x,y 2 +y 1 -l), (14) 

in terms of the inverse G of the Laplacian on the plane: 



G(x,y) 



a ixk 1 +iyk 2 



d 2 k 



4-7T 2 4 — 2 cos ki — 2 cos k 2 



The partition function of either ASM, defined as Z op = 
det A op and Z c \ = det A c i, clearly diverges. The same is 
true if one changes the boundary condition on a stretch 
of length n, by inserting n closed sites in an open bound- 
ary, or n open sites in a closed boundary. However one 
should expect that ratios be well-defined. Let us de- 
note by Z op (n) = det A op (n) the partition function of the 
ASM with open boundary condition everywhere on the 
boundary except on n adjacent sites, which are closed, 
and by Z c \{n) = det A c \(n) the partition function for the 
inverse situation. 

We consider the ratios Z op (n)/Z op and Z c \(n)jZ c \. 
They give the expectation value of the closing or opening 
of n sites on the boundary, and so should correspond in 
the scaling regime to the 2-point function of the bound- 
ary changing operator discussed above. On the basis of 
CFT, we therefore expect that in the large n limit, 



Z op (n) Z A {n) 



ZqT. 



(0(O)#n)) 



CFT 



An* , 



(16) 



for A a normalization constant. 

The explicit calculations described below confirm this 
limit, and the CFT picture behind it. 



is the determinant of a dimension n matrix. Due to the 
horizontal translation invariance of G op , the determinant 
has the Toeplitz form det (di-j), with entries a m = 5 m ,o — 
G op ((0, 1); (m, 1)) along diagonals. Using Eq. (|l3|) , a little 
of algebra gives a m as Fourier coefficients of the function 
(see Appendix A of Q for material related to the lattice 
Green function) 



a(k) = v(3 — cos k) (1 — cos k) + cos k — 1. 



(18) 



The Toeplitz determinant det (aj-j) may be computed 
by using a generalization of the Szego limit theorem due 
to Widom |||. Let cr(fc) = r(fc)(2 - 2cosfc) Q , with 
a > — i and r a single- valued, smooth, and nowhere 
vanishing nor divergent function on the unit circle, have 
the Fourier coefficients s m = ^e~ zkm a(k). Then one 
has the asymptotic value of the Toeplitz determinant 

det(s l _ J )i< lJ <n ~ E[t, a] n a2 e nt0 , n > 1. (19) 
The constant 

a G(a + l) 2 



(15) E[t, a] = e 



-(o)]- 



G(2a + 1) 



(20) 



is explicit and given in terms of the Fourier coefficients 
t m of log r and the Barnes function G . 

Comparing with (|l^), Widom's theorem applies with 

a = \ and r(fc) = \J ^ c 2 osk - J 1 ~ c 2 osk - The Fourier 
coefficients of logr(A;) are related to values of the inverse 
Laplacian on the plane ( |l5| ) 

form^O, 

1 2^[G(m- 1,0)- G(m + 1,0)] for m ^ 0. ( > 

The asymptotic ratio of the ASM partition functions 
is thus 



Z op {n) 
Z p 



E[r- I 



e 



(22) 



The exponential factor is due to the smaller free en- 
ergy that a closed site has with respect to an open 
site, a difference equal to — , as we have seen earlier. 
This term must be subtracted before comparing with 
the corresponding CFT quantity. This leaves the ex- 
pected power law with exponent 4 . The normalization is 
A = E[t; \] = e £m>i G(§) 2 - 1.18894, for which 

we have used the functional relation G(z+1) — T(z)G(z) 
and G(l) = 1 E|. 



Closed sites in open boundary 

From ([l2]), the operator A op (n) only differs from A op 
by the diagonal entries corresponding to the set I of the 
n sites which are being closed, A op (n) = A op — B, with 



It follows that 



= det (I - G op B) = det (I - G 



opJijG/ 



(17) 



Open sites in closed boundary 

In this case, the ratio of partition functions, 
Z c \(n) 



Z r 



det (I + G cX B) = det (I + G c i) i>jeI 



(23) 



poses an immediate problem: all entries in this determi- 
nant are infinite, due to the divergent integral (111]). From 
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the ASM point of view, the divergence is a manifestation 
of non local effects and is caused by the breach we make 
in the boundary, which in effect frees an infinite number 
of configurations, which were until then forbidden and 
which now become recurrent. The divergence is present 
for n = 1 and does not sharpen when n increases. 
This suggests to consider instead the ratio 



Z c \(n) 



1 



1<2 j <n i 



3 — cos fc 



Zd(l) b 

where the entries b m = 5 m ,o 
Fourier coefficients of cr'(k) — , 

v ' 2 2 y 1— cos k 

plained above, only the ratio is well-defined, but the 
numerator and denominator may be regularized by re- 
placing a 1 by 



(24) 

G cl ((0, l);(m, 1)) are the 
As ex- 



<r' a {k) = |(1 — cos k) a [yl — cos fc + V3 — cos fc] . 



(25) 



For a > — \ , all terms in (| 



are finite, but both bo and 
the determinant develop a simple pole at a — — i, so 
that the ratio is regular. 

The singularity of bo is easily evaluated, yielding bo = 
?7r(a+i/2) + nn ite. From (fl9|), the pole of the determinant 
is all contained in the factor G(2a + 1 



-1 _ r(2a+l) 
G(2a+2) 



. The coefficients t' m may be computed at a = 



1 — cos k 



3 — cos k 



IS 



(2a+l) 

where the function r'(fc) — y — ^ — , y — ^ — — 7(fey 
the inverse of the function encountered in the case "closed 
in open" , so that the coefficients t' m = —t 
those given in (^lj). Altogether one obtains 



are minus 



Z c \{n) 
ZA1) 



= 27rRes(£[T'; a]; a = e 



(26) 



By using (^o|), the prefactor is actually equal to E[t; |], 
on account ofG(|)=r(|)G(|), and produces the same 
normalization constant A for the boundary changing op- 
erator. As expected, the exponential factor now reflects 
the excess of boundary free energy of open boundary 
sites. 



We should mention that the divergent factor 



gol(l) 



has 

a finite part, related to that of G(0, 0). Its value however 
depends on the shape of the rectangle, and thus has no 
well-defined infinite volume limit. 



THE FOUR-POINT FUNCTION 

One may go further and consider the insertion of two 
stretches of closed sites in an open boundary (or vice- 
versa) , and then compare it with the 4-point function of 
0ci,op_ rpjjg con f orma i blocks of the latter are the com- 
plete elliptic integrals K(x) and A"(l — x) The ap- 
propriate solution here is (z^ = z, — Zj, and x — |— f 24 ) 



A 2 [z 12 z 34 }i 



-(1 



x 4 



K{x). 



(27) 
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FIG. 1: Ratio of the ASM partition functions for two segments 
of closed sites ([21,22] and [23,24]) in an open boundary, as 
function of the anharmonic ratio x. The solid curve is the 
CFT prediction, the other two are numerical. 



This may be directly compared with a ratio of ASM par- 
tition functions. We consider an open boundary on which 
we close two sets of boundary sites, I\ and I2, ranging 
respectively over [21,22] and [23,24]. Then the ratio 



A 2 



Z12 Z34 4 



Z p(Ii, I2) z op 

Z p{Il) Z ov {l2) 



(28) 



is expected to converge to ( |27| ) in the scaling limit. 

Exact ASM calculations are more difficult in this case, 
and we merely present in Fig.l the results of (modest) 
numerical calculations. We have taken the two sets of 
closed sites to have equal length n and to lie N sites apart 
(i.e. 221 = 243 = n, 2 32 = N). We have fixed In + N to 
200 and 300 (resp. long and short dashes), and in each 
case, we have let n run so as to make x = ( - ? Ar ) 2 vary 
between and 1. The numerical evaluation of (|28|) yields 
the dashed curves, while the solid line is the CFT result 
(P7[) . The agreement is satisfactory in the region where 
the scaling regime is best approached. 
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